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ABSTRACT 

Research is described on the fluid dynamic behavior of high performance diffusers for 
centrifugal compressors, with particular application to small gas turbine engine applications. 

Using a unique swirl generator, experiments have been carried out to define the performance and 
stall onset behavior of a modem discrete passage diffuser as a function of inlet conditions. Two 
diffusers were examined, one with 30 passages and one with 38 passages. Inlet blockage and axial 
asymmetry were varied over Mach numbers up to unity and over a range of inlet swirl angles. 
Diffuser pressure recovery and operating range were calculated using traverse measurements 
made upstream of the diffuser. It was found that the performance of the diffuser under different 
inlet conditions could be expressed to a high degree of accuracy as a single curve of non- 
dimensional static pressure recovery coefficient, based on availability-averaged inlet stagnation 
pressure, and momentum-averaged inlet flow angle. Unsteady pressure measurements showed 
that the diffuser entered rotating stall at reduced flow rates. No long wavelength stall precursor 
was determined from the measurements. 

INTRODUCTION AND PROBLEM STATEMENT 

This document represents the final report on research sponsored by the Army Research 
Organization under Grant DAAL03-90-G-0138. The research addresses the fluid dynamic 
behavior of high performance discrete passage diffusers for centrifugal compressors, with 
particular application to small gas turbine engine applications. Experimental results are presented 
for a 30-passage and a 38-passage diffuser designed at an aircraft engine company. 

Background 

The diffuser of a high performance centrifugal compressor often limits the machine’s 
useful operating range due to the onset of aerodynamic instability, in the form of rotating stall or 
surge, at mass flow rates below design. Measurement of flow conditions at the inlet to the diffuser 
is difficult due to flow field sensitivity to instrumentation. The flow field of a centrifugal 
compressor is thus typically measured at the impeller inlet and at the diffuser exit to provide 
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designers with information about the impeller-diffuser combination. These measurements, 
however, do not isolate the individual component performances. 

Defining the impeller and diffuser performance is desirable in design as well as when 
trying to optimize a tested machine by redesign. This performance split is often estimated based 
on some combination of static testing, calculated impeller and diffuser performances, and 
performance of previous designs, but little data appears to exist which shows clearly the effect of 
inlet conditions on diffusers. Testing a cascade passage diffuser in a manner in which the inlet 
conditions can be controlled is difficult due to the transonic inlet Mach numbers and high swirl 
angles characteristic of the impeller exit flow field. To address this problem, a test facility was 
constructed which enables determination of diffuser performance up to near sonic conditions 
(where the flow field is especially sensitive to instrumentation) and over a range of independently 
variable inlet flow conditions. 

The results of the investigation have been reported in depth in Refs. [1] and [2]. For this 
reason the present report will not describe details of the design of the facility, nor give an extensive 
compilation of the results of the experiments. Rather, we will present an overview of the facility 
capability, and highlight the results which appear to be of most interest. The reader is referred to 
[1] or [2] for additional material which supplements this report. 

Re s e arch O bje ct i ves 

One objective of the present study is to provide description of the causal links between 
centrifugal compressor diffuser behavior and inlet conditions. A second is to characterize the 
conditions that determine the limit of diffuser stable flow range. The variables of interest included 
inlet Mach number, inlet angle, inlet mass flow distribution, and instrumentation blockage. A 
third objective is to define the unsteady processes that are associated with the onset of instability 
(rotating stall or surge) of the diffuser flow field. 





EXPERIMENTAL APPARATUS 
Test Facility 
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The central elements of the test apparatus consists of a swirl generator which delivers an 
axisymmetric swirling transonic flow into the test section where the radial diffuser is installed. 
Reasons for the use of this method, as well as design goals, are given in Ref. [1]. The diffuser exit 
flow is dumped to a plenum and routed to a slave compressor, which can be activated to lower the 
diffuser back pressure. The overall scheme is depicted in Figure 1. The swirl generator, shown in 
Figure 2, consists of an inlet, a negative reaction impeller composed of 71 high solidity blades 
(Figure 3), a variable speed drive, and four slot rings for boundary layer control by mass injection 
or suction. Both suction and injection may be applied in an axisymmetric manner upstream and 
downstream of the impeller on either wall. The swirl generator operation envelope includes a 
maximum average exit Mach number of unity, and a range of (momentum) average inlet swirl 
angles up to 75 degrees. 

In the description below, the axial direction is taken from the impeller hub, aft wall where 
x = 0, to the shroud forward wall where x = b. The positive tangential direction is clockwise aft 
looking forward, and zero degrees at top-dead-center. “Plane 1” is the swirl generator exit radius, 
which is the traverse probe radius. “Plane 2” is the diffuser exit radius upstream of the dump to 
the plenum. “Plane 1A” is the radius defined by the diffuser inlet circumferential static taps. 
“Plane 0” is the impeller exit. All other planes are explicitly called out. 

Ra d ial Diffuser figa msiry 

There are four contiguous parts of the tested passage diffusers: the quasi-vaneless space, 
the throat, the diffuser, and the dump. Passage diffusers are defined along a centerline that is 
inclined at a large angle from a radial line. The quasi-vaneless space, which is designed to diffuse 
supersonic flow, is defined by the scalloped leading edges which span the passage width on either 
side of the passage (see Figure 4). The area converges to the throat, which fixes the maximum 
flow rate through the diffuser. The throat section is a cylinder with a length approximately equal to 
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its radius. Downstream of the throat, there is an expanding section which is a partial cone; the 
cone is cut by the two walls of diffuser which linearly diverge from the throat to the exit. The 
diffuser is thus a mix between a conical and a two dimensional diffuser, where most (over 85% 
for the 38-passage diffuser) of the area increase occurs in the tangential direction about the passage 
centerline. Table 1 lists the dimensions of the two diffusers used in the study. 


TABLE 1: DIFFUSER DIMENSIONS 


Item 

30-passage Diffuser 

38-passage Diffuser 

Inlet width 

Inlet radius 

Exit radius 

Area ratio 

Diffuser L/D 

Diffuser throat diameter 

Diffuser choking flow 

0.354 inches 

7.982 inches 

11.046 inches 

4.29 

8.75 

0.5057 inches 

2.07 pps 

0.354 inches 

7.982 inches 

11.046 inches 

4.37 

9.44 

0.4493 inches 

2.07 pps 


Instrumentation 

The swirl generator output was measured by a cylindrical, single hole probe, as shown in 
Figure 5, installed at a fixed tangential location, 0 degrees, at the diffuser inlet. The hole was 
moved axially to obtain the swirl angle and total pressure at fifteen points. At each axial station the 
probe was rotated to determine the pressure distribution on the cylinder. A second order 
polynomial was fit through the pressure distribution to determine the flow angle; fifth order and 
moving third order curve fits were used to obtain total pressure. 

One passage, for each diffuser, had a line of seventeen static pressure taps down the 
centerline, with two in the throat region, as shown in Figure 6. There were also eight extra taps in 
the quasi vaneless space, and three additional taps at the diffuser exit. Twelve static taps, six on 
each wall, were placed circumferentially around the inlet of the diffuser. One static probe was 
inserted at the same radius as the traverse probe. 

One high speed pressure transducer was placed in the plenum. Three high speed static 
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pressure transducers were installed in the vaneless space ring for the 30-passage diffuser testing. 
Five more were added for the 38-passage diffuser testing, although two of the original set were 
found to be inoperative. 

The temperature was measured near the facility inlet screen and at one location in the 
plenum. A downstream venturi was used to measure mass flow. 


DEFINITION OF PARAMETERS 


Average Values 

The inlet axial flow distribution could be controlled to achieve different levels of distortion 
velocity. The availability-averaged total pressure, as defined by Livesey and Hugh [3] and 
Filipenco [1] 



f ln(Ptj(x)) pj(x)V rl (x) 2rtr 1 dx 
JO 

fb 

Pl(x)V rl (x)27Crjdx 



was used to define the diffuser pressure recovery coefficient. 


Cpr y = 


PS 2 - Ps, 


Pt^-Ps, 


( 2 ) 


The diffuser average inlet swirl angle was found by taking the ratio of tangential and radial 
velocities mass weighted means, where 


Inl¬ 


and 




J Q p 1 (x)V rl (x)V 0 (x)2m- 1 dx 
J q p 1 (x)V rI (x)27tr 1 dx 

J Q Pi(x)V rl (x)V rI (x)27tr 1 dx 
J* Pi(x)V r] (x)27tr]dx 

vU 


(3) 


( 4 ) 


9 


This averaging process implies a momentum weighted swirl angle. 



The diffuser inlet Mach number was calculated using Pt^j 







(5) 


( 6 ) 


Profile Descriptors 

The inlet profiles of mass, momentum, energy, and swirl angle can be defined in terms of 
distortion parameters, which indicate the departure from uniformity, and skew parameters, which 
indicate the axial position about which the flow quality is equally divided. The inlet profile mass 
distortion parameter is 


X, b [(P V r)max-P<*)V r (x)] ?/ir,dx 

0,n ^iMpVrU, 

The inlet mass skew parameter is: 



b/2 


where the location of the mass split line, x* m , satisfies 


(7) 


( 8 ) 


f m p(x)V (x)2rcr,dx = f * p(x)V r (x)2Ttr 1 dx . 
Jo Jx m 


(9) 


Other parameters can be defined for momentum, kinetic energy and swirl angle, as given 
by Filipenco [1] or Johnston [2]. For example, the swirl angle distortion parameter used was the 
rms angle deviation from the momentum-averaged angle. 


Blockage Definitions 

Inlet blockage is a significant parameter for channel diffuser performance. For example. 
Sovran and Klomp [4] and Dolan and Runstadler (5] present data for two-dimensional and conical 
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diffusers which show decreasing diffuser static pressure coefficient with increasing boundary layer 
blockage at the diffuser inlet. In the present experiment, the boundary layer thickness is measured 
at the passage diffuser quasi-vaneless space inlet, not at the inlet of the channel part of the diffuser. 
In addition, for off-design high incidence cases, the boundary layer blockage associated with the 
flow around the scalloped leading edge may have a strong effect on overall diffuser performance. 

As given by Dolan and Runstadler [5], the throat blockage is defined as 

B = l _fei£tua!) . (10) 

l m ideal /throat 


The throat blockage is not directly measured by this experiment. The actual mass flow is 
measured by the venturi. The ideal mass flow can be estimated as 

Mj cos a] 


* - fT fti A i 

m ideal — \j R 




( 11 ) 


The blockage at the quasi-vaneless space inlet is based on the mass flow from Eq. (11) using the 
availability-averaged total pressure (instead of the upstream reservoir total pressure used by Dolan 
and Runstadler). In addition, it is assumed that the flow is isentropic from the inlet to the throat, 
that the throat static pressure is equal to the lower of the two measured throat static pressures, and 
that the flow is essentially uni-directional at the throat. 


TEST PLAN 

A series of constant corrected speed data sets (speed lines) were obtained for both 
diffusers. Boundary layer control was used to vary distortion parameter levels. Representative 
inlet swirl angle profiles for the no boundary layer control case and for two high skew cases are 
shown in Figure 7. The corresponding absolute Mach number profiles are also different from 
each other because the Mach number decreases as the swirl angle increases. The 30-passage 
diffuser testing had negative skew profiles, while the 38-passage diffuser testing included both 
positive and negative skews. A view of the range of inlet conditions is given in Figure 8 which 
shows the mass skew parameter as a function of mass distortion parameter for all data points, and 
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in Figure 9 (for the 38-passage diffuser) which shows the flow angle versus Mach number 
operating points for different suction and injection conditions. Each speed line consisted of a 
number of steady state operating points from choke to stall, including a near stall point. Three 
corrected speeds were examined: 2 KRPM, 4 KRPM, and 6 KRPM. Stall was detected by 
observing the output of one high speed pressure transducer, but it was also quite evident from 
decreased venturi mass flow and the plenum blow down. At the conclusion of the 38-passage 
diffuser test the traverse probe was removed and the 4 KRPM nominal and high distortion profiles 
speed lines were repeated. Unsteady data was measured without boundary layer control at several 
corrected speeds. Data was taken with and without the traverse probe being installed. 

The flow delivered by the swirl generator was axisymmetric to a high degree. The inlet 
static pressure circumferential variation was typically less than 2% of the dynamic pressure, 
pVf, with the most variation being 4%. The two independent measurements of mass flow, the 
downstream venturi mass flow and the integrated inlet profile mass flow, agreed to within 5% for 
both diffuser tests; it may be noted that 0.5 degree error in the traverse angle would change the 
integrated inlet mass flow by 3.3% for a near stall swirl angle. 

By applying suction aft of the impeller, it was possible to move the relative operating 
points of the diffuser and impeller separately. This technique was used to determine that it was the 
diffuser that initially encountered stall instability. 

RESULTS 
Steadv-State Data 

As discussed by Filipenco [1], a diffuser pressure recovery based on a dynamic pressure 
that used the availability total pressure provided a more general measure of diffuser performance 
than alternate dynamic pressure definitions (e.g. dynamic pressures based on the peak profile total 
pressure, area-averaged total pressure, or area-averaged velocity). At stall, the availability total 
pressure based diffuser recovery, Cpr v , decreased from 0.72 to 0.68 as the Mach number was 
increased from 0.15 to 0.80 under nominal conditions for the 30-passage diffuser, and to 0.67 at a 
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Mach number of 0.95. 

Pressure recovery coefficients versus the inlet swirl angle for six speed lines, three without 
distortion and three with high distortion, are shown in Figure 10 The differences in speed 
(equivalent to Mach number) and profile distortion are collapsed by using the pressure coefficient 
based on availability total pressure; alternative definitions of recovery would lead to diffuser 
performance that varied versus speed and inlet distortion. 

Pressure recovery coefficients for 20 near stall points with the 38-passage diffuser are 
shown as a function of inlet Mach number in Figure 11. (Note the expanded scales.) The diffuser 
pressure recovery coefficient decreases as the Mach number increases. No corresponding trend is 
observed for diffuser performance versus inlet mass distortion parameter, as shown in Figure 12. 
There is a weak trend of decreasing diffuser performance as the mass is skewed towards the 
impeller shroud, as given in Figure 13. In Figure 13, there are six pairs of points which comprise 
near stall points, at similar Mach number and inlet distortion parameter levels, which feature mass 
distributions that are reflected about the diffuser axial centerline, x/b = 0.5. Five of the six pairs 
indicate that when the mass distribution is skewed towards the impeller shroud, the diffuser does 
not perform as well. The one exception is where the near stall point for the negative skew case 
was not as close to stall, in terms of throttle closure, as the positive skew case. The extreme pairs, 
where suction and injection are applied, at 2 KRPM and 4 KRPM best illustrate the trend. The 
total data set of near stall points is shown in Figures 14 and 15 for both diffusers; the same general 
trends are exhibited for both. 

The question of why diffuser performance is not sensitive to inlet distortion parameter is an 
issue of interest. The relation between the (calculated) blockage at the throat and the inlet mass 
distortion parameter is definitively not a strong one, particularly at off-design, as shown in Figure 
16. We are currently examining this point; the idea is to “separate” the diffusion up to the throat 
from the diffusion downstream of the throat viewing the latter as similar to a channel diffuser. 

Removing the traverse probe had a large effect on stall onset point. Near stall, axial motion 
of the traverse cylinder caused the diffuser to stall on repeated occasions. Typically, a near stall 
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operating point was set and the machine allowed to stabilize. The test sequence would begin and 
the traverse would start to move axially with the consequence that the diffuser would enter rotating 
stall immediately after the probe moved. To determine the effect of the circumferential distortion 
created by the traverse cylinder upon diffuser stall margin, the probe was removed and 4 KRPM 
speed lines with nominal inlet and both extreme skews were measured. 

Without the traverse there is some uncertainty in determining the extent of inlet distortion 
for the skewed inlet distributions. Comparisons between the nominal inlet operating points are 
thus more meaningful but the mass flow operating range for both high skew cases increased by 
approximately 10% when the probe was not installed. 

Removing the traverse probe increased the mass flow range by 15%, with the nominal inlet 
profile. Figure 17 shows the average diffuser inlet static pressure at plane 1A versus venturi mass 
flow for the 4 KRPM throttle for three cases: 30-passage diffuser with the traverse probe installed, 
38-passage diffuser with the traverse probe installed, and 38-passage diffuser without the traverse 
probe installed. 

No direct data is available for the 30-passage diffuser performance without the traverse 
probe. Only the throttle positions and the plane 1 static pressures are available, but these do indicate 
a similar mass flow range increase. At 4 KRPM, both diffusers stalled at a plane 1 pressure ratio 
of 0.985 with the traverse installed and both stalled at a plane 1 pressure ratio of 1.0 without the 
probe installed. 

Indirect evidence is furnished from the 30-passage diffuser data taken with a bent traverse 
probe (i.e., larger blockage). The stalling flow for this configuration is estimated to be 30% higher 
than for the configuration without the probe (at a plane 1 pressure ratio of 0.95). The implication 
is that the diffuser may be more sensitive to steady circumferential distortion than to axial 
distortion. 

Based on the measured swirl generator performance, the relationship of pressure recovery 
to venturi mass flow can be derived assuming that the flow in the vaneless space is isentropic. 
Utilizing this approach, the pressure recovery of the diffuser was estimated for the nominal inlet 
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operating points taken without the traverse probe installed, as shown in Figure 18. There is a 0.04 
difference in pressure recovery, which is significant when viewed in terms of the length increase 
needed to obtain the same increase in peak recovery. A six degree conical diffuser with a peak 
pressure recovery coefficient of 0.71 would have to be 63% longer to reach the 0.75 pressure 
recovery level (Ref. [5] Figure 56, Mach number = 0.4, throat blockage = 0.03). 

Unsteady Flow Phenomena 

Unsteady measurements were also taken to define the phenomena at stall onset. 
Intermittent small amplitude, long wavelength activity was observed before stall with no traverse 
probe installed. These waves moved in the opposite sense to the impeller rotation and the rotating 
stall waves. No coherent activity was evident prior to stall when the probe was installed, as seen in 
a typical time trace shown in Figure 19. The rotating stall instability was always first observed in 
the quadrant immediately in the direction of impeller rotation “downstream” of the probe. Prior to 
stall, the static pressure variations in this region were larger than in the other three quadrants. 
Removal of the probe eliminated the probe noise at the same throttle setting, as seen in Figure 20. 
With no probe, there was coherent activity on an intermittent basis prior to stall, as seen in Figure 
21. Most of the energy was in the first and second mode, and appeared near the impeller rotation 
frequency. The inception process was short in duration, lasting about five rotor revolution times 
from first observation to essentially fully developed. 

SUMMARY AND CONCLUSIONS 

A facility has been designed and constructed for examination of diffusers for high 
performance centrifugal compressors. The facility provides a transonic swirling flow at diffuser 
inlet, with control of inlet velocity profiles. Using this facility, the fluid dynamic behavior of two 
modem discrete passage diffusers has been assessed with primary focus on the effect of inlet 
conditions. 

Definition of a pressure rise coefficient using an availability-averaged total pressure 
provided a generalized view of diffuser data for different Mach numbers and inlet distortion. The 
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mass-averaged total pressure was close to the availability-averaged total pressure, but area- 
averaged total pressure was quite different, i.e. Cpr^ = Cpr ma ss * Cpr^a • 

There appears to be little relation between blockage, or mass distortion parameter, at the 
quasi-vaneless space inlet to the throat blockage. Efforts are now ongoing to define the relation 
between inlet quantities and the throat blockage; the latter is viewed as the main factor in the 
channel diffuser performance. 

The effect of the probe on the stall onset implies that the stalls with the probe inserted are 
due to the probe wake. (The probe blockage is roughly 10% in terms of the probe cross-sectional 
area compared to the 38-passage diffuser throat area. Once the probe was removed the operational 
range increased and the stall mechanism changed. 

No periodic activity was measured at the throttle setting where the diffuser stalled when 
the probe was inserted. As the throttle was closed, small amplitude, intermittent waves travelled in 
the opposite direction to impeller rotation were detected. Rotating stall inception occurred quickly, 
within five impeller revolutions of first evidence. 
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Figure 2: Swirl generator schematic. 







Figure 3: Impeller blade shapes. 
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Figure 4; 38-passage diffuser geometry sketch. 
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Figure 5: Traverse Probe. 
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Figure 7: Variation in inlet angle versus axial position at 4 KRPM for three cases: 
O no boundary layer control, Q aft suction and forward injection, 

A forward suction and aft suction. 
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SUMMARY OF ALL PASSAGE DIFFUSER TEST DATA 



Figure 8: Summary plot of diffuser inlet conditions examined. 





MOMENTUM AVERAGED INLET SWIRL ANGLE (deg) 


25 



2K NOMINAL 
2KFW0SUC 
2K AFT SUC 
2K BOTH SUC 
2KFWDINJ 
2KAFTINJ 
2K BOTH INJ 
2kSUOINJ 
2k INJ/SUC 
4K NOMINAL 
4KFWDSUC 
4K AFT SUC 
4K BOTH SUC 
4KFWD INJ 
4K AFT INJ 
4K BOTH INJ 
4k SUC'INJ 
4k INJ/SUC 
6K NOMINAL 
6K BOTH SUC 


Figure 9: Flow conditions examined (38-passage diffuser): inlet swirl angle and inlet Mach 

number combinations. Legend indicates speed and inlet control, where K = 1 KRPM, 
inj = injection, sue = suction, FWD = shroud side (x/b = 1), and AFT = hub side 
(x/b = 0). 
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Figure 10: Diffuser overall, availability-averaged pressure recovery, Cpy, as a function of 
momentum-averaged inlet swirl angle (30-passage diffiiser). 
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Figure 11: Near stall operating points, Cpy versus Ml for 38-passage diffuser. Legend indicates 
speed and inlet control, where k = 1 KRPM, inj = injection, sue = suction, FWD = 
shroud side (x/b = 1), and AFT = hub side (x/b = 0). 
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Figure 12: Pressure recovery coefficient at near stall points versus inlet mass flow distortion 
parameter. Legend indicates speed and inlet control, where K = 1 KRPM, inj = 
injection, sue = suction, FWD « shroud side (x/b = 1), and AFT = hub side (x/b = 0). 
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Figure 13: Pressure recovery coefficient at near stall points versus mass flow skew parameter. 

Legend indicates speed and inlet control, where K = 1 KRPM, inj = injection, sue = 
suction, FWD = shroud side (x/b = 1), and AFT = hub side (x/b = 0). 
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Figure 14: Comparison of diffuser near stall Cfty versus inlet Mach number for the 30-passage 
and the 38-passage diffusers. Nominal inlet conditions are indicated by circles, 
suction or injection is indicated by squares, and 38-passage data is indicated by shaded 
symbols. 
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Figure 15: Comparison of near stall Cpy versus inlet swirl angle for the 30-passage and the 38- 
passage diffusers. Nomina] inlet conditions are indicated by circles, suction or 
injection is indicated by squares, and 38-passage data is indicated by shaded symbols. 
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Figure 17: Diffuser mass flow range with and without traverse probe. 
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THROAT BLOCKAGE VERSUS INLET BLOCKAGE 



Figure 16: A comparison of the estimated throat blockage and the inlet blockage. 
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ESTIMATED CP FOR DIFFUSER WITH AND WITHOUT PROBE 



Figure 18: Estimated 4 KRPM speed line pressure rise coefficient for nominal inlet conditions, 
with and without traverse probe. 






4000 RPM w/PROBE 



high speed pressure data throttle to stall. Elapsed time was one second. 
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w/o PROBE 
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Nominal inlet, no traverse probe, 4 KRPM high speed pressure datp throttle to stall. Elapsed time was one second. 
Normalized time, x, is t^i/27tri. 






















